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Abstract 

We have studied the electron transport through a quantum dot cou- 
pled to three leads in the presence of external microwave fields supplied 
to different parts of the considered mesoscopic system. Additionally, we 
introduced a possible nonresonant tunneling channels between leads. 
The quantum dot charge and currents were determined in terms of 
the appropriate evolution operator matrix elements and under the wide 
band limit the analytical formulas for time-averaged currents and differ- 
ential conductance were obtained. We have also examined the response 
of the considered system on the rectangular-pulse modulation imposed 
on different quantum dot-leads barriers as well as the time-dependence 
of currents flowing in response to suddenly removed (or included) con- 
nection of a quantum dot with one of the leads. 



1 Introduction 

The electron transport via resonant tunnelling in mesoscopic systems has been 
the subject of extensive theoretical research due to recent development in fab- 
rication of small electronic devices and their potential applications. Some 
interest has been focused on the transport properties of a quantum dot (QD) 
under the influence of external time- dependent fields. New effects have been 
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observed and theoretically described, e.g., photon-assisted tunnelling through 
small quantum dots, [H |2] , photon-electron pumps |2| and others. In most 
theoretical investigations a QD placed between two leads was considered (e.g., 
Refs. 121111130 El El El) and the current flowing through a QD under periodic 
modulation of the QD electronic structure or periodic (nonperiodic) modu- 
lation of the tunnelling barriers and electron energy levels of both (left and 
right) leads was calculated. 

One of the important problems of the mesoscopic physics is the interference 
of the charge carries. This interference appears when two (or more) transmis- 
sion channels for electron tunnelling exist. Such possibility exists, e.g. in the 
electron transport through a QD embedded in a ring in the Aharonov-Bohm 
geometry and much theoretical interest has been paid to description of the 
phase coherence in this and related systems, e.g. Refs. pUl ITTl [T2\ . Another 
experimental situation in which the interference may occur can be realized 
in the scanning tunnelling microscope (STM). The recent experimental and 
theoretical studies with a low-temperature STM of a single atom deposited on 
a metallic surface showed the asymmetric Fano resonances in the tunnelling 
spectra, e.g. Refs. [311211111 ■ In the STM measurements tip probes the trans- 
mission of electrons either through the adsorbed atom or directly from the 
surface. The transport of electrons through both channels leads to an asym- 
metric shape of the conductance curves which is typical behavior for the Fano 
resonance resulting from constructive and destructive interference processes. 
The quantum interference can be also observed in the mesoscopic system with 
multiple energy levels [15^ . A model which incorporates a weak direct nonreso- 
nant transmission through a QD as well as the resonant tunneling channel was 
also considered in Ref. JH] in the context of the large value of the transmission 
phase found in the experiment for the Kondo regime of a QD [T7j . 

A number of works has been devoted to the problem of the electron trans- 
port in the multiterminal QD systems and here we mention only a few of them. 
In Ref. PH] the conductance of the N-lead system was considered showing that 
the Kondo resonance at equilibrium is split into N-1 peaks. In Ref. ^H] an 
explicit form for the transmission coefficient in the electron transport through 
a QD connected with three leads was found. The electron transport and shot 
noise in a multi-terminal coupled QD system in which each lead was disturbed 
by classical microwave fields were studied in Ref. Multiterminal QD sys- 
tems or magnetic junctions were also intensively investigated in context of the 
spin-dependent transport, e.g. Refs. [211 123 • A three-terminal QD system 
was studied in Ref. [221 121] to measure of the nonequilibrium QD density of 
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Figure 1: Schematic picture of the QD coupled with three leads. It can serve 
as a possible STM experimental setup when the left lead (L) corresponds to 
the STM tip. 

states (splitting of the Kondo resonance peak). The cross correlations of the 
currents and the differential conductance of the QD coupled with three leads 
described by the infinite-U Anderson Hamiltonian were considered in Ref. |25j . 
The general formulation of the time- dependent spin-polarized transport in a 
system consisting of the resonant tunnelling structure coupled with several 
magnetic terminals was considered by Zhao et al. [22] and as an application 
of this formalism the electron transport in a system with two terminals under 
an ac external field was investigated. 

In all papers mentioned above and in the studies relating to the electron 
transport through a QD with the additional (bridge) nonresonant transmission 
channels, the time-dependent external fields were not applied and the consid- 
ered systems were driven out of the equilibrium only by means of a dc voltage 
bias (see, however, Ref. [22]) 

In this paper, we address the issue of a QD coupled with three leads with 
additional, nonresonant coupling between leads driven out of equilibrium by 
means of a dc voltage bias and time- dependent external fields. The QD is 
connected with three metal leads and one of these leads, say the left (L) lead, 
is coupled with the two remaining leads, say the first (Ri) and the second 
{R2) right leads. The possible STM experimental setup corresponding to our 
model system is presented in Fig. 1. In literature, different theoretical ap- 
proaches have been developed to treat the time-dependent electron transport 
in the mesoscopic systems. The most popular nonequilibrium Green's func- 
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tion method depends on the two time arguments and for time-dependent prob- 
lems it is a rather difficuh task to calculate them without any approximations 
(e.g. beyond the wide-band limit). Therefore, in our treatment of the time- 
dependent problems we use the evolution operator which, as a rule, essentially 
depends on one time argument (e.g. Refs. |2Z1 12H1 123)- Such an approach is 
especially well suited for the problem with time-dependent coupling between 
the QD and leads. 

The outline of the paper is as follows. In Sec. II we start with the model and 
method for the derivation of the QD charge and currents. In Sec. Ill we present 
the results for the time-averaged currents and their derivatives with respect to 
the QD energy level position (or equivalently, with respect to the gate voltage) 
obtained for different time-dependence of the parameters characterizing the 
considered system. We consider also the transient current characteristics in 
the case of the rectangular-pulse modulations imposed on the QD-lead barriers. 
The last section presents the conclusions and in the Appendix we give the short 
derivation of the evolution operator matrix elements needed in the QD charge 
and current calculations. 

2 Model and formalism 

We consider a QD coupled through the tunneling barriers V^^ {i = 1,2,3) 
to three metal leads. One of these leads, say the left lead (L) is coupled 
additionally with the remaining two leads, say the first and second right leads 
(-Ri, -R2) by the tunneling barriers Vj:^ In the following we will denote 

the wave vectors associated with the left lead by the letter k and the wave 
vectors corresponding to the first and second right leads by the letters q and 
r, respectively. The chemical potentials /ij of the three metal leads may not be 
equal, and their difference is not necessarily small. We write the Hamiltonian 
of the considered system in the form H = Hq + V, where 

Ho= £p{t)a^aj^ + edit)ajad, (1) 

p=k,r,q 

V = iVp^d{t)at ad + h.c) + 

p=k,f,q 

+ E(l^^,-(t)ata, + h.c.) + $:(l^^.(t)ata,-.+ h.c.) (2) 

k,r k,q 
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The operators ap-(ai), ad{a^) are the annihilation (creation) operators of the 
electrons in the corresponding leads and the dot, respectively. For simplicity 
the dot is characterised only by a single level Ed and the intra-dot electron- 
electron Coulomb interaction is neglected. The neglect of Coulomb interaction 
is quite reasonable in some systems and, as we are going to concentrate on the 
investigations of the influence of the third lead (in comaprison with the QD- 
two leads system) and the additional tunneling channels between the leads 
on the time- dependent transport, then in the first step ignoring the Coulomb 
interaction should be justified. We consider our mesoscopic system in the pres- 
ence of external microwave (MW) fields which are applied to the dot and three 
leads. In most theoretical treatments of photon-assisted electron tunneling it 
is assumed that the driving field equals the applied external field. However, 
the situation is more complicated and the internal potential can be different 
from the applied potential [H]. One of the consequences will be, e.g. the 
asymmetry between the corresponding left and right sidebands [201 120] • The 
main feature of the time-dependent transport remains, however, unchanged 
and in our treatment as usual we assume that in the adiabatic approxima- 
tion the energy levels of the leads and QD are driven with the frequency u 
and the amplitudes Aj {i = L,Ri,R2), and read £j:-(t) = + Ajcoscut, 
£d{t) = £(1 + cosa;i:, respectively. 

The time-evolution of the QD charge and the current flowing in the system 
can be described in terms of the time-evolution operator U (t, to) given by the 
equation of motion (in the interaction representation): 

z^^^ = vm{t,to), (3) 

with V{t) = Uo{t, to) V{t)Ufl{t, to) and Uo{t, to) = Texp (^i / dti Ho{ti)^ where 

T denotes the time ordering and the units such that h = 1 have been chosen. 
Here we have assumed that the interactions between the QD and leads, as well 
as between the left and right leads are switched on in the distant past to- 
The QD charge is calculated according to the formula (cf. Refs. fI7\ 

ndit)=ndito)\Uddit,to)f+ nj}ito)\Udpit,to)f- (4) 

p=k,r,q 

Here Uddit, to) and Udp(t, to) denote the matrix elements of the evolution oper- 
ator calculated within the basis functions containing the single-particle func- 
tions \k), \q), |f) and \d) corresponding to three leads and QD, respectively. 
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ridito) and np{to) represent the initial filling of the corresponding single-particle 

states. 

The tunneling current flowing, e.g. from the left lead, Jl(^), can be obtained 
using the time derivative of the total number of electrons in the left lead, 
jL{t) — —edniit) / dt (in the following we put e = 1), where 

Mt) = EM^) = EM^o)it/fe,(t,to)r+ 

k k 

+ E%,(^o)i%,(t,to)r+ 

k,ki 

+ E^f(*o)|C/fe,-(^,io)|' + E^Kio)IM^.^o)r. (5) 

k,q k,r 

In the following only the matrix elements of the evolution operator present 
in Eqs. (4) and (5) are required and they can be obtained solving the corre- 
sponding sets of coupled differential equations constructed according to Eq. 
(3) with Vab{t) written as follows: 

Vab{t) = Vab{t)exp (^i{ea - eb){t - to) + ^ ^'' (sina;^ - sina;to)^ , (6) 

where a and b correspond to |^), q), \r) or \d), respectively. 

As the example, the matrix element Udd{t,tQ) required for the calculation 
of the first term of the QD charge, Eq. 4, can be obtained solving the following 
set of coupled differential equations 

« E ydp{t)Upd{t,to) , (7) 

p=k,q,r 

^Vuamddit, to) - ^ E %(^) Up<iit^ to) , (8) 

p=q,f 

tVpdmdd{t,to)-^J2%iiit)UkditM, P = q,r. (9) 

k 

The total number of coupled equations in this case is equal to (3A^ -|- 1), 

— * 

being the number of discrete wave vectors k, q and f taken to perform the cor- 
responding summation over the wave vectors. Usually, the number N equal 
to 100-200 is sufficient to achieve the desired accuracy of the calculations. We 
have solved numerically this and other similar sets of the coupled differential 



dUdd{t,to) 
dt 

dt 

dUfdjt, to) 
dt 
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equations needed in calculations of all matrix elements of the evolution oper- 
ator present in Eqs. 4,5. We have used this method for the special case of 
time-dependent couplings of the QD with leads and the couplings of the left 
lead with two right leads. The set of Eqs. 7-9 in the case of vanishing over-dot 
couplings between the left and right leads is greatly simplified and gives, e.g. 
for Udd{t,to): 



dUddit,to) 
dt 



' dhK{t,h)Udd{t,to) , (10) 
to 



where 



K{t,h) = ydpitWpditi) = 

p=k,q,r 

J2 \Vi\\i{t)ui{ti)Di{t - ti) exp{iedit - h)) 

i=L,Ri,R2 

X exp(i(Arf — Aj) (sin cut — sinti;ti)/co') (11) 

and Dj(t — ti) is the Fourier transform of the z-th lead density of states and 
Vdi{t) = ViUi{t). Similar simplifications occur in the calculations of other ma- 
trix elements of U{t, to) required in the formulas for nd{t) and ni{t). However, 
for the nonvanishing couplings V^^ and V^- (over-dot bridge between the left 
and right leads) one has to solve the starting set of Eqs. 7-9. Much more 
analytical calculations can be done in the case of constant values of the tun- 
neling matrix elements present in our model. In this case the general equation 
satisfied by Udd{t, to) is derived in the Appendix and under the wide band limit 
(WBL) approximation, e.g. [21 IH El this equation takes the simple form 

^^^^ = -C^if/..(Mo), (12) 

here C*! = (f - T^lf^) T, x = ttVrl/D, D being the bandwidth of the lead 
energy band (-D/ji = = = D) and V = 2ttV'^/D. In the Appendix we 
give the derivations of all functions needed for calculation of the QD charge 
and currents. We assumed the simplified assumption that all tunneling matrix 
elements are independent of the wave vectors. The interactions between the 
QD and leads are assumed to be equal between themselves and denoted by 
V and the interactions between the left and two rights leads (i.e. V^^ and 
V^^) corresponding to the over-dot tunneling channels are also equal one with 
another and denoted by Vrl. 
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It is easy to show that the first term of the general formula for the QD 
charge, Eq. 4, together with the solution of Eq. (12), Udd(t, to) = exp{—Ci(t — 
to)), tends to zero for t — to — ^ cxd as ReCi = ^ — i^2x^ ^ ^- '^^^ T^ext terms 
of the QD charge formula can be calculated using the functions U^^it.to), 
Udqit,to) and Udrit^to), Eqs. 43,47, being the solutions of the corresponding 
differential equations, Eqs. 42,46. Finally, the time-averaged QD charge is 
given by 

{Mt)) = E / def,{e){\A,{e,t)\^) , (13) 



=L,Ri,R2 



where 



= (1 + 4x^)7(1 + 2a;')2r/27r, (14) 

(15) 



UL = l-L -t- )/{L -t- ) i /ZTT , 

aji^ = qr^ = (1 + x^)/(l + 2x^)^1/271 . 
t 

Ai{e, t) = -i j rftiexp {-i{ed - e){t - ti) - z(Ad - A^) 

to 

(sin^t-sin^ti)/^)) exp ( ^^^-^^-^(t-ti)j , (16) 

where (...) denotes the time-averaging and fi{s) denotes the Fermi function of 
the i-th {i = L,i?i,i?2) lead. Noticing, that lm{Ai{e,t)) = - j^(|Ai(e, t)|2) 
(cf. 0), the expression for the time-averaged QD charge can be written as: 

(1 j_ 47.2 p 

+ ^;^^E//«,(^)(.4*(M))<fe). (17) 

In order to calculate the current ^^(t) the functions Ui^^(t,to), Ui^_^i^^{t,tQ), 
U^^it^to) and Uj^^(t,to) are required and they are given in the Appendix in 
Eqs. 37, 44, 49. After lengthly but straightforward calculations we obtain for 
the time averaged current leaving the left lead the following formula: 



i=Ri,R2 



2x2 , 



7r(l + 2x2)2 

+ g(^I fL{e){AL{e,t))de- I f,{e){A{e,t))de 
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where 



G 



37r(l + 2x2)3 



(Qx{1 - 2x^) + i(l - 13a;^ + Ax 



(19) 



A,, - A, 



a; 



e — Ed — (^k 



2Tx . 3r/2 ^ 
1 + 2x2 ^ * 1 + 2x2 



(20) 



and Jfe(y) denotes the Bessel function. The corresponding formula for the 
time-averaged current {jRi{t)) leaving i?j-lead, i = 1,2, cannot be written in 
such symmetrical form as in Eq. 18, because the i?j-lead is coupled with L-lead 
only. For {jRi{t)) we have: 



+ 



2x' 



7r(l + x2)2 

r 



37r(l + 2,x2)3 V^^' / fR.i^)i^Rii^^t))de-G, J h{e){Ar^{e,t))de 



G, I fR.{e){AR.{e,t))de 



where Gi = 6x(l - 2^2) + i(l - 13^2 + 4x^), G2 = 3x 



(21) 

(-2 + 5x2 ^ ^4)^ 

G3 = 12x=^ + i(l + 8x2 _ 5^4-) ■ ^ 1^2) for i = 2(1). Note, that the in- 
tegrals present in the formula for the QD charge and currents, Eqs. 17,18,21, 
can be easily performed analytically and final algebraic expressions can be 
obtained. Especially simple and transparent form can be given for the con- 
ductance -^{jjit)), i,j = L,Ri,R2. For example, 



reads as: 



_d_ 



7r(l + 2x2)2 



j:4 

k 



A, -A, 



^ /r2(l - 13x2 + 4x4) ^Yx{\ - 2x2) N 



7r(l + 2x2)^ 



where 



/ 2rx \^ / 



7r(l + 2x2)s 



3r/2 



(22) 



+ 2x2 



(23) 



2rx 

1 + 2x2 



(24) 
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Analyzing Eq. 22 one can find the origin of tlie asymmetric line shapes in differ- 
ential conductance resulting from the interference of resonant and nonresonant 
tunneling paths. For the case of Vlr = we observe the Lorentzian resonances 
localized at Ed = /^l ±uk. The amplitudes of these resonances are determined 
by the fc-th order Bessel functions calculated for the argument (A^ — Ai)/uj. 
For the case of nonvanishing Vlr, the resulting curve is a superposition of 
the Lorentzian-like resonances and asymmetric parts weighed by the factors 
1^(1 -13a;2 + 4x^)/(l + 2x2)^ and ATx\l-2x'^)/'K{l + 2x'^f , respectively. The 
Lorentzian-like resonance is centered at Sd = /iL±Ci;A; + 2rx/(l + 2x^) with the 
peak width at half maximum (FWHM) equal to SF/ (l+2a:^) and the maximum 

value equal to -^^-^^^^p-Jl ( ^''^^^ ) • The asymmetric part of the differential 
conductance corresponding to the k—th sideband is also centered in the same 
point with the distance between its maximum and minimum equal to and 

the absolute values of these extrema are equal to ^ li+2x^)^ '^k ( ^''~^^ )- For 
comparison, in the case of the QD coupled with two leads the corresponding 
Lorentzian-like part of the differential conductance corresponding to the fc— th 
sideband is centered at Sd = /iL ± cu/c + Fx/(1 + x^), with FWHM equal to 

2F/(H-x^) and the maximum value equal to ^ ^(1+2x^)2^ "^fc ( ^''^^^ )- Knowing 
the explicit expressions for the currents one can check the following relations 
between different elements of the conductance matrix —ed{jn(t)) /dfim, e.g. 
[IHllSni- Current conservation implies J2nd{jn{t)) /dfim = 0, 
On the other hand, J2m9{jn{t)) /dfim = only for - A^^ = - Ar^ = 
Ad — A/jj. For other relations between the amphtudes A^ and A^, A^^, A^^ 
we have: 

j:d{jRAt))/df^k = E^O^.W)/^/^^ = -lY.d{jLit))/df^, (25) 

k k k 

for Ad - Ar, =Ad- Ar, ^Ad- Al, and 

Y.d{jRAt))/dfik ^T.d{jR.{t))/dfi, ^Y.d{jL{t))/dfik (26) 
for Ad - Ar, ^ Ad - A^, ^ Ad - A^. 
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3 Results and discussion 



We consider the QD coupled with three, say the left, the first right and second 
right metal leads with the additional over-dot (bridge) couplings between the 
left and right leads. We present the results for the time-averaged currents and 
derivatives of the average current with respect to the QD energy level in the 
presence of external microwave fields which are applied to the dot and three 
leads, respectively. The time-dependent currents are also calculated in the 
case when the periodic rectangular-pulse external fields are applied to each 
QD-lead barrier. 

We have taken for Vlr the values comparable with ^ and estimated 

^kad (assuming its /c-independence, Vg^^ = Va = V) using the relation Ta ~ 
SttIVqP/Z^o, where Da is the a-lead bandwidth and Da = 100 Ta (F^ = Tr = 
T, Dl = Dfi = D was assumed). In our calculations we assumed e = 1, all 
energies are given in F units, time in ^/F units, the current, its derivatives 
and frequency are given in eT/h, e^V/h and V/h units, respectively. 

Firstly, we consider the time-averaged currents in the presence of the time- 
varying (harmonic case) external fields. Here we give the explicit formula for 
the averaged current, (jl(^)), performing the corresponding integrals in the 
general formula given in Eq. (18), and compare it with the current {jl {t)) 
flowing in the system of the QD coupled with two leads only: 



{Ut)) 




(1 + 2x2)2 



2^2 



F 1 - 13x2 + Ax^ 
37r (1 + 2x2)3 




Fx(l -2x2) 



TT (1 + 2x2)3 



X E[4^(^^)ln(..)-J.^(^^)ln(,0]), 



(27) 



where hi = (//j -ea-uk^ T^S^ ) / ( 2(1+2^2) ) ' 9i = {l^i -Sd-u;k + 



l+2a;2 



2Tx 





and i — L,Ri,R2, whereas for {j^ (t)) we have: 



2x2 



2 (/^L - f^R) + 



F 1 - 6x2 + x^ 
2^ (1 + 2x2)3 



7r(l + x2) 
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Figure 2: The averaged current {ihit)) against Sd in the system of a QD 
coupled with three (sohd hne) or two (broken hne) leads. The microwave field 
is applied to the QD with = 6 and = 5, /iL = 0.2, iir^ = = 0, F = 1, 
Vrl = 0. 



X 



+ 



X 



E 

k 



Jk { — ~ ~^ arctan(/i^^'') — ^ — arctan(/i^'*^ 



TT (1 + 

E 



k 



{2i 



where hf' = (/i^ -ed-ujk+ / ( (h^) and ^ff ^ = (/i^ - - cj/c + 

(^i+x^)2 , i = L, R. Note, that {jhit)) consists of the two terms and each term 
is similar in its structure to the current fiowing in the QD-two leads (QD- 
2L) system, {jl {t)). However, due to the interference of the charge carriers 
propagating along the different ways the arguments of the arctangens and log- 
arytmic functions are different and the individual terms in (jl(^)) and (j^ (t)) 
are weighted by different x-dependent factors. 

In Figure |21 we compare the averaged values of the current flowing from 
the left lead in the systems in which the QD is coupled with three or two 
leads, the solid and broken curves, respectively. The external microwave fleld 
is applied only to the QD and dc bias between the left and right leads is small 
in comparison with u, and F. The coupling Vlr is assumed to be zero. In 
such the sidebands on the current curve are clearly visible. The width 
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Ed 

Figure 3: The averaged current {jiit)) against Ed in the system of a QD 
coupled with three (thick hues) or two (thin hues) leads. The solid, broken 
and dotted curves correspond to x = 0, 0.28 and 0.7, respectively. = 8, 
Aifi = = Al = 0, fiL = -fJ'Ri = -/iija = 0.1, u; = 5. 

of the corresponding peaks is smaller for the case of the QD coupled with two 
leads. Analyzing the expressions for {jiit)) and {jL'\t)) one can obtain (for 
/iL << r and /iRi = /ii?2 = 0) the subsequent peaks in {ibif)) (as functions of 




with the FWHM equal 3r. For 



comparison, in the case of the QD coupled with two leads the corresponding 

peaks are described by the functions -^Jl ('^) A^l ^1 + (^^^y^) ) with the 
FWHM equal 21. 

Next, in Fig. El we analyze the influence of the over-dot coupling between 
the left and right leads. Again, as in Fig. we consider an external microwave 
field apphed only to the QD with = 8 and Hl = —fJ'Ri = —f^R2- The thin 
(thick) lines correspond to the QD coupled with two (three) leads. We show 
{jiit)) for three values of the inter- leads coupling strength represented by the 
parameter x = ttVlr/D. For x = we have well defined sideband peaks as 
previously shown in Fig. |21 Next, we show the results for x = y^ ^^-vTsf ^ g 2g 
and X = \/2/2 ^ 0.7. Note, that the expression for (jiit)) (and (jl^''(t))) 
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consists of three terms. The first term depends on the difference jj^i — Hi and 
does not depend on the QD energy level Sd- The second term corresponds to 
the Loretzian-type contribution to a given sideband (it disappears for x ~ 0.28) 
and the last term corresponds to an asymmetric contribution (it disappears 
for X ~ 0.7), respectively. This last term influences the sideband shape only 
for nonvanishing over-dot coupling between leads and is the most prominent 
sign of the interference effects. For x = V2/2 this term disappears and the 
resulting sidebands have a form of symmetric dips due to the negative value of 
the coefficient ^£~/_j^^^+3^^ in Eq. 27. In this case the nonresonant tunneling 
channels modify the dips center position and its FWHM in comparison with 
the position and FWHM of sidebands presented for x = 0. Note, that for a QD 
coupled with two leads, the corresponding sidebands (the thin dotted line) are 
not fully symmetric curves as in this case the last term of Eq. 28 disappears 
for X = 1 and not for x = 0.7 (c.f. ref. |23). For x ~ 0.28 the corresponding 
sidebands are fully asymmetric curves as in this case the second term in Eq. 
27 (which introduces asymmetry) disappears. Again, for a QD coupled with 
two leads the corresponding sidebands (the thin broken curve) are described 
by not fully asymmetric curves as the second term of Eq. 28, which gives a 

symmetric contribution to sidebands, disappears for x = — 2^/2 ^ 0.38 
and not for x = 0.28. 

In the next Fig. 0] we show the currents flowing in the QD-three leads 
(QD-3L) system for a; = and x = 0.28. For vanishing values of the over-dot 
coupling the current (jiit)) is characterized by a sequence of the symmetric 
peaks, but the current (jflal^)) ^ superposition of the asymmetric structures 
placed in the points where the symmetric sidebands occur on {jiit)) curve. 
Analyzing the current (jRiit)) according to Eq. 29 we have for the parameters 
in Fig. m 

(j/?2(^)) = ^J2Jk (—) [arctan(/;,ijj - arctan(/ii) 

+ arctan(/;,ij2) — arctan(/;,/jJ] , (29) 

where hi = {fii — Cd — ujk) One can see that each sideband is the sum of 
the peak (two first terms in Eq. 29) and the dip (the last two terms in Eq. 
29) resulting in the asymmetric structure shown in Fig. |3] 

To learn more about the influence of the third electrode and additional 
over-dot coupling between leads we present in Fig. El the currents {jhif)) for 
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Figure 4: The averaged currents {ihif)), {jRiif)) and (ii?2(^)) (solid, broken 
and dotted lines, respectively) against Sd in the system of a QD coupled with 
three leads. The upper (lower) panel corresponds to x = (x = 0.28). = 8 
and other parameters as in Fig. El The values of {jR2(t)) for x = (upper 
panel) have been multiplied by a factor of 20 for illustrating purposes. 

X = and x = 0.28 and {jRzit)) for x = 0. For vanishing coupling between 
leads (x = 0) the current {jLit)) exhibits a well known sideband structure for 
T < u and for small frequencies, a; < F, the two broad maxima at Ed = ±Ad 
are present. At the same time, the current 0i?2(^)) exhibits the asymmetric 
structures centered on the (e^, u) plane at the points where photon sidebands 
occur on the {jhif)) curves. These asymmetric structures on the {jR2(t)) curve 
exist also at u < T. On the other hand, similar structure of the (jiit)) (the 
lower panel of Fig. is obtained for x = 0.28, i.e. we observe a number of 
asymmetric resonances separated by the photon energy for a; > F. Notice the 
similarity of both pictures, i.e. {iiit)) calculated for x = 0.28 and {jR2it)) for 
X = 0, respectively. Note, however, the different scale for these currents. 

In Fig. El we present {jiit)) for different over-dot coupling assuming a 
strong asymmetry of the applied microwave field (ac potential is applied only 
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Figure 5: The averaged currents against Ed and u. The upper (middle) panel 
shows {ihit)) {{jR2(t))) for X = and the lower panel gives (jL(t)) for x = 0.28. 
Ad = 8, Ar^ = Ar^ = Al = 0, = -/iiJi = -/iiJa = 0-2- 

to the i?i-lead in the QD-3L system). The additional i?2-lead is characterized 
by the chemical potential fiR^ = (a^l + /^Ki)/2 = 0. For better demonstration 
of the influence of the ac potential on {jL{t)) we moved down each curve by 
the constant value A = 2x^{2^i — fiR-^ — fiR^)/{l + 2x'^Y, see Eq. 27. For 
X = we observe a shoulder on the left side of the main peak and a small 
negative current for the positive values of Ed- This picture is similar to the 
known results (for x = 0) obtained experimentally and theoretically in the 
QD-2L systems, e.g. Refs. |31 EEj- With the increasing over-dot coupling 
Vlr the height of the main resonant peak decreases and disappears at all for 
X = 0.28 but the enhancement of the current for negative values of Ed (for 
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Figure 6: The averaged current (jiit)) against Ed for different values of 
the over-dot couphng strength between the left and right leads. The bro- 
ken, dotted, thin solid and thick solid lines correspond to x = 0,0.14,0.28 
and 0.7, respectively. A = 2x^(2^1 — /i/j^ — + 2x^Y ^-Ri = 3, 

= Ajj2 = Ad = 0, = -/i/j, = 0.2, = 0, u; = 5. 

~ ~5) is more and more expressive. At the same time, for all values of x 
we observe a negative current for the small positive Ed with a greater absolute 
value for stronger coupling between leads. For greater values of x the shape 
of the {jiit)) curve is changed dramatically and for x = 0.7, {ihif)) becomes 
nearly reversed in comparison with that calculated for x = 0. 

In Fig. [3 we analyze the influence of the third lead (here named as R2) on 
the current when the external microwave field is applied to this lead 

and to the QD with A^ = 2/S.ji^. For comparison, we add in the upper panel 
the results for {ihit)) obtained for the case when this additional third lead is 
not irradiated by the microwave field. In this before, see e.g. Fig. 

121 we observe typical sidebands on the current curves (the difference between 
the lead chemical potentials is small in comparison with the amplitude A^). 
However, after including the third lead irradiated by the external microwave 
field the dependence of the current {jLif)) on the gate voltage (or equivalently 
on the QD energy level position) is quite different - compare the thin or thick 
solid lines of the upper and lower panels. For smaller values of A^^ and , the 
averaged current is very similar to the corresponding current Jl obtained 
by applying the external microwave field only to one lead (see the broken 
line in Fig. IHl These curves are, however, related between themselves by 
a relation Jii^Sd) — jL{—£d)- Now we can observe a small negative current 
at small negative values of Ed and some enhancement of the current on the 
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Figure 7: The averaged current {ibif)) against Sd for the case of the microwave 
field apphed only to the QD with Ad = 3 (thin line) or = 6 (thick line) - 
upper panel. The lower panel corresponds to the case of the microwave field 
applied to the QD and R2 lead with Ad = ^Ar^ and Ar^ = 3 or 6 (thin or thick 
lines, respectively). The broken lines show the results when the microwave field 
applied to the QD and R2 lead are out of phase (with the phase difference of 

T^), f^L = -fJ'Ri = 0.2, =0, UJ = 5. 

right side of the main peak. Similarly, the significant differences between the 
corresponding currents are observed also for greater values of the amplitudes 
Ad and Ar^ (compare the thick solid lines in the upper and lower panels of 
Fig. Ej). Note, that very similar behavior of the current (jiit)) as the function 
of the gate voltage is observed if we compare the case when the microwave field 
is applied only to one lead and the case when the microwave field is applied 
simultaneously to the QD and i?2— lead but with the phase difference of vr - 
compare the thin broken curve in Fig. [7|with the broken curve in Fig. IHl 

Next, we consider the QD-2L system for different values of the chemical 
potential of the right lead (yU^j = —0.1,4 and 10) and show the results of in- 
cluding the next lead (i?i— lead) to the system with /i/j^ = —0.1. In Fig. |Hlwe 
present the results obtained for (jl(^)) in the case when the microwave field is 
applied only to the QD. For the case fi^ = —f^Ri = —f^R2 = 0-1 ^^nd x = the 
sidebands are very clearly visible and the corresponding peaks are lower and 
broader for the QD-3L system as we discussed before (the upper panel). For 
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Figure 8: The averaged current {jL{t)) against Sd in the system of a QD 
coupled with three (thick hues) or two (thin hues) leads. The upper (lower) 
panel corresponds to a; = (x = 0.28). In a QD-three leads system fiL = 0.1, 
/iRi = —0.1 and fiji2 = —0.1,4 and 10 (solid, broken and dotted curves, 
respectively). In a QD coupled with two leads /il = 0.1, fiR = —0.1,4 and 
10 (solid, broken and dotted curves, respectively). The values of solid curves 
have been multiplied by a factor of 10 for the illustrating purposes, u = 5, 
Aa = 8, Ar, = Ar, = = 0. 

the non-zero over-dot coupling between the left and right leads the sideband 
peaks get asymmetric forms and for x = 0.28 they become fully asymmet- 
ric (the lower panel). The current for the QD-2L system is also composed 
of a number of asymmetric components although now these forms are fully 
asymmetric for x = 0.38 as we know from the earlier discussion. For greater 
values of /irj, the corresponding current flowing in the QD-2L system achieves 
(for vanishing, as well as for nonzero coupling between leads) greater negative 
values and its dependence on the QD energy level position is well-marked in 
comparison to the results characterizing the QD-2L system. 

In order to emphasize the influence of the additional lead on the currents 
we show in Fig. Elthe current {jL{t)) for QD-2L and QD-3L systems calculated 
for the parameters for which the corresponding curves are relatively simple. 



19 



-10 -5 5 10 
Ed 



Figure 9: The averaged current {jhit)) against Sd for the QD-2L systems: 
Hl = —fJ'R = 0.1 - thin sohd hues and = 0.1, = 4 - thick sohd hues 
and the QD-3L system: /xl = — = 0.1, fiR^ = 4 - broken hues. The 
upper (fower) part corresponds to x = (a; = 0.28) and = 5, = 1, 
A/j^ = A/J2 = Al = 0. The dotted hues correspond to the sum of the currents 
flowing in two QD-2L systems. 

We assumed the small amplitude of the QD energy level oscillations, A^ = 1, 
and a; = 5, for which (for x = 0) only the central peak corresponding to elastic 
tunneling is visible on the {jiif)) curves. We show {jiit)) obtained for two 
different QD-2L systems which can be viewed as components of more com- 
plicated QD-3L system. We observe that due to the interference effects, the 
current {jiit)) flowing in the QD-3L system is not simply a sum of currents 
(dotted lines) flowing in the corresponding QD-2L systems. The difference 
between this sum and the current corresponding to the QD-3L system is rela- 
tively large and exists independently of the coupling between leads. 

In order to present more information about the differences between the 
electron transport in the QD-3L and QD-2L systems we display in Fig. 
the derivative d{jL{t)) /dsd as a function of Ed and u. The lowest panel corre- 
sponds to the QD-3L system and two other panels correspond to the QD-2L 
systems. These QD-2L systems are characterized by such parameters that 
combined together give us the considered QD-3L system. One can see, that 
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Figure 10: The averaged current derivatives d{jL{t)) / ded against Ed and uj for 
a; = 0. We compare the results obtained for the QD coupled with two leads 
- the upper (middle) panel - for = 5, /xr = —8, A/, = 8, = 4, Aij = 
(/iL = 5, /iij = 0, Ai = 8, Arf = 4, Ar = 2) with the results obtained for the 
QD coupled with three leads - the lower panel - Aj^ = 8, A^ = 4, A^^^ = 2, 
Ar^ = 0, /ii = 5, /iij, = 0, /iR2 = -8, u; = 5. 

the considered characteristics of the electron transport in the QD-3L system 
are not simply the algebraic sum of the corresponding curves of both QD-2L 
systems. In all three cases shown in Fig^J the position of the corresponding 
minima and maxima (along the e^— axis) can be identified with the values of 
the leads chemical potentials. However, the corresponding structures are less 
clear in the case of the QD-3L system in comparison with those for the QD-2L 
models. 

In the next step of our investigations of the electron transport in the QD- 
3L systems we consider the time-dependent currents flowing in response to 
the time-dependent barriers between the QD and leads or in response to sud- 
denly removed (or inserted) connection of the QD with one of the leads. As 
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Figure 11: The time-dependent current flowing in the system of a QD coupled 
with the three leads: L, i?i and i?2- The L-lead is coupled with the QD only 
- the left panels and with the QD and two other leads, Vlr^ = Vlr2 = 4 - 
the right panels. The couplings between the QD and Ri, R2 leads are changed 
periodically. The upper, middle and lower panels correspond to /i/, = 3, and 
-3, respectively. —fJ.R^ = ixr^ = 3,6^ = 0. The thin, thick and broken curves 
correspond to JljJri and Jr^ currents, respectively. 



a first we assume a rectangular-pulse modulation applied to the QD-Ri lead 
and QD-i?2 lead barriers. Wc assume that these modulations are with a phase 
difference of tt. In the first (second) half-cycle, V^r^ = (V^/jj = 0) ^-iid the 
QD is coupled only to the R2 lead {Ri lead). In addition, the QD is coupled 
to the next, say L— lead, with a constant value V^l- In the following we con- 
sider the time-dependent currents jiit), Jri(^) and jR^it) for the three specific 
conditions: fJ,L — /J'Ri, IJ'L — (iJ'Ri + IJ^R'z)/'^ ^^nd hl — A* Hi. In addition, we 
assume ^ir,^ = — /Xri = 3, Ed = Uri and take for the period of the consid- 
ered barrier modulation T = 5. In a such case we integrate numerically the 
corresponding set of the differential equations for the matrix elements of the 
evolution operator and in the next step calculate the currents according to 
the formula ja{t) — —edna{t)/dt, where a — L,Ri or R2 and na{t) is given in 
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Eq. 5 (or similar to it). We checked that the QD charge hardly depends on 
the additional V^/j couplings. Although the QD charge is almost insensitive to 
the additional over-dot couplings the currents demonstrate such dependence. 
Especially visible are the differences for the case when the chemical potential 
Hl of the third electrode L lies between the chemical potentials of two other 
leads, see Fig. ^2 B,E. For other values of fi^, the influence of the over-dot 
tunneling channels for the parameters considered here is smaller. Note, that 
after abrupt changing of the coupling between the QD and Ri, or R2 leads 
the currents jL,jRi and j/jj rapidly change too, and after a short time reach 
the steady values. The QD coupled with three leads could be considered as 
the three-state system. We observe that for some values of the lead chemical 
potentials the currents change their values periodically e.g. from zero to the 
positive value (see Jr^ in Fig- IHB), from zero to the negative value (see Jr-^ in 
Fig. Illb) or from the negative to the positive value (see Jl in Fig. Illbl.The 
additional couplings between leads modify the values of the currents but the 
qualitative picture remains the same. Note, that in the first moment after 
abrupt changing of the coupling between the QD and Ri or R2 lead we have 
jLit) + j_Ri(t) + jR2(t) 7^ as in this case dnd(t)/dt ^ (not shown here). 
After some delay time the QD charge stabilizes around its equilibrium value, 
the currents tend to constant values and their sum is equal to zero. 

In the last step we consider the response of the currents to the abrupt 
inclusion into the QD-2L system of the third electrode (in our case, L-lead). 
The results are presented in Fig. together with the schematic view of 
the subsequent tunnel connections between the QD and the three leads. We 
show the time-dependent currents jL(t), jRiit) and jR2it) corresponding to 
the three different ways of inclusion of the L-lead. We assumed the chemical 
potentials fiR-^ = —fiR2, I^l = (/^i?i + /^i?2)/2 and Ed = yUija- Consider the 
current jr^ flowing from the i?2-lead characterized by the highest chemical 
potential hr^. Before adding to the system of the L-lead the current jr^ has 
a constant value and flows from i?2-lead through the QD energy level to Ri- 
lead. When the L-lead is included into the system (the tunneling coupling Vlr^ 
changes abruptly at t = 8 from zero to nonzero value) the current Jr^ is almost 
unchanged - its value decreases slightly without any transients at short times 
after the time t = 8. Next, we consider the case when the L-lead is abruptly 
connected simultaneously with the QD and i?i-lead. Now the current jr^ 
decreases significantly during the short time after the moment of inclusion of 
L-lead and settles to its constant value. Note, that jr^ decreases despite the 
additional charge transfer channel between R2 and L leads (through the QD 
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Figure 12: The time-dependent currents flowing from tlie L, Ri and i?2 leads 
in the system shown in Fig. ^ The thin solid, thick solid and broken curves 
correspond to different couplings of the L lead to the other elements of the 
system: {Vlr, = 4, Via = Vlr^ = 0), {Vlr, = 4, 1/^ = 4, Vlr^ = 0) and 
{Vlri = Vid = Vlr^ = 4), respectively. 



energy level). The destructive interference appears in this case as we have 
two transmission channels for tunneling electrons between R2 and Ri leads. 
However, the constructive interference is visible if we consider the next case 
when the L-lead is abruptly coupled with the QD-2L system assuming nonzero 
values of Vl/j^, Vlr^ and V^d- Now, we have one additional charge transfer 
channel (from i?2-lead) in comparison with the former case. The current Jr^ 
rapidly increases with some fluctuations and after the time ^ h/T decreases 
to the constant value greater than in QD-2L system. Similar analysis can 
be made considering the currents ji and Jr^ although the transient current 
changes are more visible now at short times after the abrupt inclusion of 
additional electron tunneling channels. The above discussion concerns the 
specific values of the lead chemical potentials and the position of the QD 
energy level. Nevertheless, the similar qualitative conclusions can be made 
also for other values characterizing the considered system. 
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4 Conclusions 



We have studied the time-dependent tunnehng transport through the QD cou- 
pled with three metal leads using the evolution operator technique. The time- 
dependent QD charge and currents were determined in terms of the appro- 
priate evolution operator matrix elements. Applying the wide band limit to 
the intcgrodifferential equations satisfied by the evolution operator matrix el- 
ements we were able to give the analytical expressions for the time-averaged 
currents and differential conductance. We considered the external harmonic 
microwave fields applied to different parts of the considered system, as well 
as, the rectangular-pulse modulation imposed on different QD-leads barriers. 
In addition, we have studied also the time dependence of the currents due 
to abrupt inclusion into the QD-two leads system of the third electrode. We 
have considered also the effect of the additional couplings between leads (we 
coupled one of the leads with the other two leads) on the conductance and 
current fiowing in the system. Our main results can be summarized as follows: 

- For the vanishing nonresonant timneling path, Vir = 0, and for the param- 
eters for which the photon-assisted sidebands are clearly visible on the {jhif)) 
curve the subsequent sideband peaks have the Lorentz-type form with the 
FWHM equal SF in comparison with 2F for the QD-two leads system. For the 
increasing value of Vlr the form of the sidebands transforms from the Lorentz- 
type to the fully asymmetric form for x — \J (13 — -\/153)/8 ~ 0.28. For greater 
X the form of the sidebands changes and gains again the Lorentz-type shape 
for X = v^/2. For the QD coupled with two leads the corresponding values of 

X are equal to y/s — a/52/3 ~ 0.38 and 1, respectively. 

- For the vanishing Vlr the differential conductance curve, e.g. d{jL{t)) /dfi^, 
possesses the sidebands of the Lorentz-type localized at //^ = ± uk. For 

Vlr 7^ these sidebands are described by the superposition of two parts, the 
Lorentz-type and the asymmetric one centered at = /xl ± ujk + 2r,T/(l + 
2x^) weighed by the factors '-^^^j^g^Ji (^) and J| (^) , 

respectively. For the QD-two leads system the corresponding k—th side- 
bands are centered at = A^l i + ^x/{l + x^) and their symmetric 
and asymmetric parts are weighted by the factors ^-^^^^jr-^ Jk ( ^'^~^'^ ) and 
?;g#J|(^), respectively. 

- The symmetry properties of the sidebands corresponding to the current fiow- 
ing from the given lead depend on the position of the chemical potential of 
this lead in comparison with the chemical potentials of the other two leads. 
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Taking, for example Vlr = and /irj localized in the middle between hl 
and II one can observe on the {jhif)) curve the sidebands of nearly regular 
(Lorentz-type) forms while the sidebands on the (j/jal^)) curve have asymmet- 
ric structures. However, in the presence of the nonresonant tunneling path the 
sidebands on the {jiit)) curve change their form and for x ~ 0.28 they have 
a fully asymmetric shape. On the other hand, the sidebands on the (jRal^)) 
curve have for x ~ 0.28 a nearly Lorentz-like shape. 

- Especially large interference effects can be observed if we compare the cur- 
rent flowing in the QD-3L system with the sum of currents flowing in the two 
QD-2L systems which can be viewed as the components of the considered more 
complicated QD-3L system. The difference between them is relatively large 
independently of the over-dot coupling between leads (see Fig. ^ 

- In the case of strong asymmetry of the apphed external field {Ar-^ 7^ 0, = 
Ad = Ar^ = 0) we observe for a; = a shoulder on the left side of the main 
resonant peak on the (jiit)) curve vs. the gate voltage. With the increasing 
over-dot coupling between L and Ri, i?2 "leads the main resonant peak disap- 
pears and transforms in a dip for strong coupling Vlr. At the same time a 
shoulder of this curve increases with the increasing Vlr. 

- Let us consider the time dependence of currents flowing in response to the 
time-dependent barriers between the QD and two leads (we assume a constant 
coupling of the QD with the third lead). For the assumed rectangular-pulse 
modulation applied to the QD-i?i lead and QD-i?2 lead barriers one can con- 
sider the QD-3L system as a three-state one. For example, the currents jL{t), 
jRiit) and jR2it) change their values periodically between zero and positive, 
positive and zero and negative and positive values, respectively (see Fig. IIB). 
The additional couplings between the leads introduce only small quantitative 
changes. 
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Appendix 

In this section we present the derivations of the general equations satis- 
fied by the required functions Udd{t,to), U^j:{t,to), Ud,f/q{t,to) and Uj:^{t,to), 
^ki ^2^^' ^o)' ^kr/q(^' ^o) nccdcd for the calculations of the QD charge nd{t) and 
the currents flowing in the considered system. In the next step, using the WBL 
approximation we simplify these equations and give the analytical solutions 
for them. Let us begin from the derivation of the integro-differential equation 
satisfied by Udd{t,to). Writing down the formal solution of Eq. (11) 

t 

Uq/fd{t,to) = -i j dtiV^/pd{ti)Udd{ti,to) 

to 

t 

- ^ J dhY.%/rkiti)UMitiM (30) 

to k 

and inserting them to the formal solution of Eq. (10) 

t t 
UkditM = -^J dhV^^{h)UM{ti,to) -ijdhY. Vkpih)Upd{ti,to) , (31) 

to to P=?^^ 

one can obtain after straightforward calculations the function U^^{t,to) ex- 
pressed in the terms of Udd{t,to), only 

t t ti 

Ukdit^to) = {-i) j dti%^{ti)Udd{tuta) + {-i? j J dhdhRj:a{tut2)Udd{t2,to) 

to to to 

+ EH)''//--- / dti...dt2j Rj^^,^{t,M) Ru^^hM) ■ ■ ■ 

i— 2 to to to ki,k2,---,kj-\ 

tj)Udd{tj, to) 

+ EH)'''^'/ /■■Jdt,...dhj+,^J2^RkkMuh)R^,^^,^{h,u)... 

to to to ki,k2,...,kj 

Vkjdi'^2j+l)Uddit2j+l,to) , (32) 

where 

Rij(h,t2) = E %(tl)VV, (t2) + E Vir{tl)Vr^j{t2) (33) 
q r 
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We remember, that the wave vectors k, q and f correspond to the left lead 
and the first and second right leads, respectively. Inserting into Eq. 9 the 
expressions for the functions U^d{t,to), Ufd{t,to) and Uj^^{t,to), Eqs. (17, 18), 
we can write the integro-differential equation for Udd{t,to) in the form 



dt 



- j dti 



to 
oo 



i 

' E ydkitWkdit^) + Rdd{t, to) Udditi, to) 

k 

+ /■■■ / dt,...dt,,^ j: jRdkSt^ti)Rk.k,it2,ts) . 

to to to ki,k2,---,kj 

^kj-lkj (^2j-2, t2j-l)Vj;.^{t2j) 

+ ^dki (*)-^fcife2(^l' ■ ■ ■ Rk.d{t2j-l,t2j)]Udd{t2j, to) 

00 * *1 *2i 

+ T.i-if' j j ■■■jdt,...dt2,+,'£[R,uStMRuS^2M)... 

to to to k\k2---kj 
^kj-lkj ('^2i-2, t2j-l) Rkjdi^2j,t2j+l) 

+ E Kifci W -^fcife2(^l' ^2) ■ ■ ■ -Rfe^,fe^+i(^2j-l, t2j)V^.^^^{t2j+l)]Udd{t2j+l: 
kj+i 

This rather untractable general equation can be greatly simplified Tising the 
WBL approximation. Assuming Vdf = Vdq = V^^ = V, = V^ = Vrl, the 
multidimensional time integrations and summations over the wave vectors can 
be performed giving in result 

dUdd{t,to) „ . s . s 
— = -Gi Udd{t, to) , (35) 

with the solution 

Uddit,to) = cxp{-Ci{t-to)). (36) 

Here Ci = |r - ^^^^^^ and x = tvVrl/D, D being the bandwidth of the lead 
energy band {Dji-^ = Dr^ = = D). The function Uj^^{t,to) given in Eq. 
(19) can be reduced within the WBL approximation to the form 

t 

Ujidit, to) = -C2 j dt,V,;,{t,)Udd{tu to) , (37) 

to 
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where 



i + 2x 



In order to calculate U^j:{t, to) we write down accordingly with Eq. (3) the 
corresponding set of coupled differential equations for the functions U^^{t,to), 
^ifcifo^^' ^o); U^{t,to) and f/-?^(t,to)- The subsequent steps of the calculations 
are similar to those performed in the derivation of Eq. (34). Inserting the 
formal solutions for Uj:j^{t,to), U^^f,^^{t,to) 



to 



dti 



p=qr 



u 



to 



ki 



(39) 
(40) 



into the differential equation satisfied by U^^{t^ to) one obtains the derivative 
dU^j^{t,to)/dt expressed in terms of U^j^{t,to) and U^_^j:^{t,to). On the other 
hand, on the basis of Eqs. (33, 34) the function U^_^^^{t, to) can be represented 
in the form containing only U^j^{t,tQ). Finally, one obtains 

t 



dU^t^to) 
dt 



dk 



dti 



to 



E^c..l%.-(^l) +E^c.n(^)VV,.-(^l) 



ri 



t ti 



+ {-i)^ J J dtidt2 

to to 

t ti t-z 

+ {-if j j j dtidt2dt. 

to to to 



E^.A;.W%,i(^i)n-.^^(^2) + E--- 

ki,qi fei,ri 



E ^^ci.l(^)%,-,(^l)%.,(^2)V-,,-(t3) + 



E ■■■+ E •••+ E ••• 

ri,kiq2 qikif2 r\k\r2 



+ ... 



+ 



i-i)^ jdt, 

to 



E VrkS^njii) + E VaqAmUt^) 



91 



+ E^^nW^nd(^l) 



UaSuto) 
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t tl 



+ 



J J dtidt2 



to to 



ki,qi 

ki,fi 91 fel rifci 



In the WBL approximation this equation reduces to the form 

dUgkit^ to) 



dt 



-C2V,^{t) + CsU,^{t,to) 



with the solution 



Udk{t,to) = -C2 1 dhV^^ih) ,ep{-C3{t - h)) 

to 



(42) 



(43) 



where C3 = r(4ia; - 3)/2(l + 2x'^). Taking into account Eqs. (33, 34) the 
function Uj^^j^^{t,tQ) can be written in terms of U^j^^{t,to) 



to 

t 



1 + 2x2 



/ 

to 



dti e 



(44) 



In order to calculate Ud^t, to) we first write down the coupled set of equa- 
tions (on the basis of Eq. (3)) satisfied by the functions Udq{t,tQ), U^^^^(t,to), 
Uff^{t, to) and Uj:^{t, to). Inserting the formal solutions for the functions U^^g-^ {t, to), 
Ufqit^to) and Uj^^{t,to) into the differential equation for the function Udq(t,to) 
one obtains the derivative dUdq{t, to)/dt expressed in terms of Udq{t, to), U^J^t, to), 
U^i^^{t,to) and Up^{t,to)- Inserting again into this equation the formal solu- 
tions for the functions U^^^^{t,to), Uf^{t,to) and Uj^^{t,to), and repeating this 
process again and again, one obtains: 



^^^^ = -K.,w + (-)^/^tiEv:.fe^iW%,(^i) 

to fel 



+ 



t tl 

(-i)^y J dtidt2 

to to 



flfel 
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t tl t2 

+ {-if J J J dtidt2dh 

to to to 



kiqik2 



felflfe2 



+ ... 



+ J dh E ydp{t)Vj;d{ti)Ud^{tuto) 

to p=ki,qi,ri 

t tl 

qi,ki 



+ 



{-if J J dtidt2 

to to 



+ E--- + E--- + E--- 

fifel kiqi fcifi 



Udq{t2,to) + ... 



(45) 



In the WBL approximation this equation is reduced to the following form: 

dUdq{t,to) 



dt 



^-C^Vd^{t)+C^Udq{t,U), 



where C4 = (1 — x)/(l + 2x^) and has the solution as follows: 

Udq{tM = -Ca f dhVdg{h)exp{-C;{t - tl)) . 

Jto 



(46) 



(47) 



The function Ucif{t,to) is identical with Udq{t,to). 

To calculate, e.g. the current jiit), we still need the functions Uy{t,tQ) 
and Uj^^(t,to). The function Uj:^{t,to) can be obtained solving the set of the 
coupled differential equations for the functions Ud^it,to), U^q-^{t,to), Uf^^i^tjto) 
and U^Jyt^to). Writing down the formal solution for Ui^^(t,to) and inserting 
into it, in the first step, the formal solutions for U^^(^^{t^tQ) and Uf^{t,to) and, 
in the second step, the formal solutions for Ud^{t, to) and U^J^t, to), and so on, 
one obtains 

t 

Ukqit.to) ^-ij dt,V^,{t,)Ud,{t,,to) 

E%i(^l)%<^(*2) +E%iW^n<i(^2) Ud,it,,to) 



to 



+ 



(-i)^ J dtidt2 

to 



11 



n 



31 



t h i2 

+ {-if J J J dtidtidts 

to to to 

t 



qi,ki 



fiki 



Udqih, to) 



t ti t2 



+ ... + i-i) I di2%(ii, to) + i-i)^ J J J dtidt2d3 

to to to to 



E%.(^i)%.^.(^2)%,(i3) + E--- 



Under the WBL approximation this equation reduces to the form: 



(48) 



Ukq^it, to) = -C2 1 dt,V^,(t,)Ud,it,, to) - ^^-^ I dt,%{t{) , (49) 

to to 

where Udq{t,to) is given in Eq. (43). The function U^p{t,to) has the identical 
form as U^J^t^to). 
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